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Abstract — A customary solution to reduce the energy con- 
sumption of wireless communication devices is to periodically 
put the radio into low-power sleep mode. A relevant problem is 
to schedule the wake-up of nodes in such a way as to ensure 
proper coordination among devices, respecting delay constraints 
while still saving energy. In this paper, we introduce a simple 
algebraic characterisation of the problem of periodic wake-up 
scheduling under both energy consumption and delay constraints. 
We demonstrate that the general problem of wake-up times 
coordination is equivalent to integer factorization and discuss 
the implications on the design of efficient scheduling algorithms. 
We then propose simple polynomial time heuristic algorithms 
that can be implemented in a distributed fashion and present a 
message complexity of the order of the number of links in the 
network. 

Numerical results are provided in order to assess the per- 
formance of the proposed techniques when applied to wireless 
sensor networks. 

Index Terms — Wireless Networks, Energy Saving, Wake-up 
Scheduling, Chinese Remainder Theorem 

I. Introduction 

In wireless networks with battery-operated devices, energy 
saving mechanisms are of paramount importance in order to 
maximize network lifetime. Measurements have shown that 
the power consumption associated to the reception of packets 
is of the same order of magnitude as of that involved in 
the packets transmission H], Q, lO. Even worse, due to 
the power drainage operated by the radio frequency (RF) 
amplifier, the energy expenditure for passive communication 
operations - e.g., overhearing the channel for collision avoid- 
ance mechanisms - has a very high energy cost as well. 
As a consequence, power saving mechanisms typically force 
low-power consumption modes for the RF interface, which 
in turn inhibits the communications capabilities (transmission, 
reception and channel sensing) of wireless devices for a certain 
fraction of time. The typical statement to this respect is that 
nodes are put into sleep mode. 

A customary solution to increase network lifetime is there- 
fore to periodically put nodes into sleep mode. This corre- 
sponds to introducing a duty cycle, intertwining sleep and 
active periods at the device level. Such approach is already 
in use in a variety of wireless technologies, including power- 
save modes and beaconing techniques in IEEE802.il H and 
wake-up cycling methods in wireless sensor networks 0. 

The introduction of wake-up scheduling has been identified 
in the past as a major source of performance degradation in 
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wireless networks. Nodes going in sleep mode increase indeed 
the latency associated to the delivery of messages. This has 
been subject to a number of research studies, in particular 
showing that duty cycles may lead to a data forwarding 
interruption problem in wireless sensor networks [0, [0. 

The problem of scheduling wake-up times may appear 
at first sight a standard optimization problem: one tries to 
maximize the wake-up periods while satisfying a given delay 
bound 17], m in order to minimize the associated energy 
consumption. However, in this trade-off there exists a further 
element of complexity since the scheduling of activity periods 
requires coordination at the network level. Indeed, nodes in 
sleep mode are not able to receive incoming messages, so 
that it is necessary that both transmitting and receiving nodes 
become active at the same time. This problem is often referred 
to as scheduling of rendezvous points among neighboring 
nodes. While it can be regarded as a synchronization issue, 
it is inherently different from clock synchronization ||9], ifTol . 
where one typically aims at synchronizing devices' clocks 
network-wide using minimal local message exchange. 

A recent research line has addressed the design of quorum- 
based protocols for the power saving duty cycles ID, ifTTI . 
lfT2l . lfT3l . IIT4II . In general, wake-up schedules should be 
robust with respect to loss of synchronization due to clock 
skew and should account for possibly different energy/delay 
constrains on various nodes (due to, e.g., different types 
of battery installed in different devices etc.). Furthermore, 
protocols for determining the wake-up schedules should be 
able to adapt to changes in network topology, due to, e.g., 
node failures or addition of new nodes to the network. In the 
aforementioned works, the target is the definition of algorithms 
and protocols able to ensure that two neighboring nodes 
are able to communicate (i.e., both are active at the same 
time) at least once within a given (bounded) time-frame. A 
different approach to reduce energy consumption is adaptive 
listening [HI, ifTSll . Q, whereby nodes enter a low -power mode 
if no activity is detected on the channel for a given time 
interval. However, in this approach, no guarantees can be given 
in terms of delay and/or energy constraints. 

In this work we take a different perspective and aim 
to answer to the following questions; Assume that energy 
and delay constraints have been assigned network-wide: do 
feasible duty-cycle schedules exist?If so, what is the related 
energy-delay trade-off? Furthermore, do algorithms exist that 
achieve network-wide distributed scheduling of duty cycles, 
and scale in the number of required messages? 

Our main contributions can be summarized as follows: 



• a simple algebraic characterization of the wake-up 
scheduling problem; 

• the analysis of the complexity of the wake-up scheduling 
problem: it is proved equivalent to that of integer factor- 
ization; 

• heuristics and algorithms that solve the wake-up schedul- 
ing problem with message complexity linear in the num- 
ber of links in the network; 

• fully distributed protocols implementing the aforemen- 
tioned algorithms and test them in a system-level simu- 
lator to evaluate their performance. 

The remainder of the paper is organized as follows. Sec. |ll] 
reports on related works in the literature. Sec. HID introduces 
the system model and formally define the problem of wake- 
up scheduling problem. Then, we analyze in Sec. |IV] the 
complexity of the aforementioned problem by relating it to 
integer factorization while Sec. FV] introduces algorithms solv- 
ing the problem with low messages overhead. In Sec. IVII we 
evaluate the performance of the proposed protocols, obtained 
by means of system-level simulations. Finally, Sec. IVIII is 
devoted to final a conclusions and pointers to promising 
research directions. 

II. Related Work 

In wireless networks, periodic wake-up is a convenient mean 
to avoid idle listening to the channel and to prolong nodes' 
Ufetime 15], ID, lH, lfT6l . Typically, the duty cycle (defined as 
the fraction of time a node is in the active state) is set to values 
of the order of some percent ||2], llT6l . This is feasible when 
the data rate is low, i.e., on the order few packets per minute, 
as in low-rate wireless sensor networks. In order to schedule 
periodic wake-ups, a simple solution, adopted by BMAC and 
other WSN protocols, is to employ preamble sampling ||5], 
O. Basically, the preamble is set long enough to guarantee 
node discovery at each wake-up. 

In order to avoid inefficiencies, it has been proposed to 
match the RF interface activation, i.e., the sleep period, to 
the traffic pattern. This is the case of SMAC or TMAC 
protocols 121, iflTl . lfT6l . In particular SMAC synchronizes the 
duty cycle by coordinating neighbouring sensors to the same 
slot to reduce the time spent listening to the channel. Adaptive 
change of the duty cycle is supported in TMAC ifTTl . whereas 
lfT6l matches duty cycles to the dynamics of the routing 
protocol. Tuning such protocols network-wide is typically not 
considered; scalability issues do indeed arise, due to the need 
of a potentially large amount of signaling messages. This is the 
core issue that we address in this paper: we design distributed 
algorithms for wake-up scheduling that present linear message 
complexity in terms of the network size and can operate 
local adaptation in case of nodes joining/leaving in a fully 
distributed fashion. 

In ad hoc networks literature, centralized planning of wake- 
up schedules was addressed in several works. Solutions that 
are close to the ones proposed in this paper can be found in , 
HD, mi, m, mi. In H the target is the design of power 
saving protocols for WiFi networks, using quorums design 



and leveraging the power-save mode provided by the IEEE 
802.11 standard. In [fTTl methods for designing asynchronous 
and heterogeneous wake-up schedules are presented, where 
the heterogeneity may come from either different applica- 
tions running or by different node-level features. Along the 
same line, lfT2l and ifTsl provide solutions for creating wake- 
up schedules whereby the ability of neighbouring nodes to 
communicate (i.e., to be active in the same time interval) 
within a wake-up cycle is ensured by combinatorial means. 
The complexity of quorum design limits the applicability of 
such combinatorial technique to regular or simple topologies; 
further insight into the design of quorums is found in fl9\ . 

Some other works address the issue of network performance 
in the presence of duty cycles fS), l?). The goal of such 
works is to maximize the throughput and minimize the delay, 
while meeting a given average power consumption constraint. 
It is shown that the problem is, in general, NP-complete; an 
analysis is provided for regular topologies (line, grid and tree 
networks in Q, tree and ring networks in ||8]). The distributed 
heuristic provided in Q does not provide guarantees in terms 
of delay bounds. Centralized approaches only are considered 
in Id. 

In II20II the authors address the joint design of routing 
and wake-up scheduling protocols. The idea is that routing 
should be aware of the wake-up schedules and vice versa. A 
distributed implementation (based on the use of distributed 
versions of Bellman-Ford algorithm) is also possible. The 
need to coordinate routing and wake-up scheduling may create 
issues related to the amount of signaling messages to be 
exchanged in rapidly changing topologies. 

Main contributions 

In this work we expose the equivalence of the wake-up 
scheduling problem to integer factorization. We observe that 
this follows naturally as soon as one imposes constraints on 
both delay and energy consumptions. Using basic algebraic 
tools 1211 we provide a framework that relies on a very 
Umited set of input parameters, providing means to map energy 
consumption constraints and end-to-end delay requirements to 
wake-up scheduling. 

We further provide fully distributed heuristic algorithms 
that solve the wake-up scheduling problems with a message 
complexity 0{M) messages, where M is the number of 
links in the network. The proposed schemes are validated 
numerically using system-level simulations. 

III. System Model 

We consider a multi-hop wireless network represented as a 
graph G{V, E) where G is a set of N vertices corresponding 
to nodes j = 1, . . . , A^, and is a set of A/ edges. Edge ij 
connects vertices j and i if the two nodes are within mutual 
communication range. We assume a discrete time model and 
denote by To the time slot duration. 

In the following we will assume Tq = 1 for notation's sake. 
We consider a common time axis, and denote the origin by 0. 



TABLE I: Notation used throughout the paper 



Symbol 




V,N 


set of nodes, number of nodes 


E,M 


set of links, number of links 


di 


degree of node i 




neighbors of node i 


H 


j-th neighbor of node i 


TV, 


time slot duration 


71; 


wake-up period at node i 


Oi 


start epoch of wake-up cycle at node i 


aij 


phase of node i towards node j 


Ai 


set of active slots for node i within a peiiod 


(r,s) 


n-th rendezvous time for nodes i and j 


greatest common divider of r, s g Z 


[r, s] 


least common multiplier of r, s S Z 



Nodes do not need to be aware of such common time axis: 
we will use it only for presentation purposes (see Sec JIV-AV 

Ng{i) is the set of neighbors of node i; we further assume 
that nodes are synchronized on time slot boundaries. 

We now define a periodic wake-up schedule. To illustrate 
the concept and notation, we refer to the simple example 
reported in Fig. [T] Two neighboring nodes, namely node 1 
and node 2, wake-up with period rii = 5 and n2 — 3, 
respectively. This means that they wake every ni (respectively: 
n2) slots. However, their wake-up cycles can start from a 
different "origin". In the example, node 1 period starts at slot 
1 and node 2 period at slot 2, respectively. In particular, let 
ai be the local time origin for node i: in the example, the two 
nodes wake up at time with respect to their local time, but 
at time 1 and 2 with respect to the common time. 

We want to characterize the conditions under which for 
every pair of nodes in the network, a wake-up schedule 
provides a set of times when nodes are both awake and can 
communicate; those times are called also rendezvous times. 

Nodes may have multiple active slots within a wake-up 
period. We denote by Ai the set of such active slots for node i. 
In other words, node i follows an activation cycle with period 
rii, whereby the active slots are characterized by the set Ai. 
We call an element of Ai a phase. 

Formally, we define a periodic wake-up schedule for node 
i as a function S{i) = {Ai, rii), where {Ai, rii) is a pair given 
by the wake-up epoch set Ai C {0, 1, . . . , — 1} and the 
wake-up period 7i.i G Z. Under a wake-up schedule, node i 
will wake up at times 

Ai + riiZ, := {a + k ■ rii, k £ 'L,a £ Ai} 

We say that node i performs a periodic wake-up with period rii 
and epochs in Ai. With a slight abuse of terminology, we will 
call a wake-up schedule the set of functions S = {Si}i=i^,,,_M\ 
also, for ease of reading, we will denote G Ai the phase 
that node i uses to communicate with node j. 

Denote X = {.t^} the set of the meeting times (or: 
rendezvous points) for nodes i and j, i.e., the set of time 
slots when node i and j are both active: X = {Ai + n^Z) n 
{Aj +njZ). As it will be clear in the following, only two 
cases are possible: X = $ ot \X\ = oo, i.e., either there are 
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Fig. 1: Wake-up schedule example: qi — 1, ni = 5, a2 — 2 and 
n2 = 3; nodes meet every 15 slots starting at slot 11. 

no meeting times or there is an infinite number thereof. Quest 
for periodicity is the reason for such dichotomy. 

We consider first energy constraints. Assuming that the 
energy expenditure is proportional to the fraction of time a 
node is in the active state, one wants to prolong the lifetime 
of node i by letting the duty cycle \Ai \/ni satisfy constraints 
of the form: 

n, > L, ■ \A,\, i = l,...,N (1) 

We consider the possibility of nodes having different energy 
constraints (i.e., different values for Li). This may reflect, e.g., 
heterogeneity in the battery capacity. 

Further, depending on the specific applications running on 
top of the network, specific requirements for the end-to-end 
delay traffic may have to be satisfied. For instance, the delay 
bound may represent the maximum latency that is allowed for 
a certain sensor reading. 

In the remainder of the paper, we assume that constraints of 
such a type are given in terms of the maximum delay between 
two subsequent rendezvous between node i and any of its 
neighbors j. We denote such constraint by L^^; the constraint 
can then be written as: 

xl+^-x^^<U„\/jeNg{t). (2) 

It is important to notice here local constraints Li and Ui 
are input to the problem. In App. |E] we show how end-to-end 
delay constraints can be turned into constraints of the form 
(|2]i, while accounting at the same time for the presence of 
energy constraints as in ([T). 

Definition 1 A wake-up schedule S is said to be feasible if 
for every node i, the number of rendezvous points with node 
j, where j € Ng{i) is infinite and if it respects the energy 
constraints ([7]). 

Definition 2 A wake-up schedule is said to be tight ;/ it is 
feasible and if it further respects the constraints Q . 

Definition 3 A wake-up schedule is said to be strictly feasible 
(respectively: tight) if it is feasible (respectively: tight) and 
\Ai\ ^lforalli^l,...,N. 

The wake-up scheduling problem can be formalized as: 

Problem 1 (Strong WAKE-UP) Given graph G = {V, E) 
and constraints (|2| and (Q, determine a (strictly) tight wake- 
up schedule. 



A trivial necessary condition is that Ui > L + 1 = 1 + 
max{Ljj, . . . , ij„}, yijk S E, which we assume holds in 
the remainder of the paper 

If a wake-up schedule is feasible but not tight, some of the 
constraints (|2]l are not satisfied. This origins a weaker problem 
that is 

Problem 2 (Weak WAKE-UP) Given graph G = {V,E) 
and constraints (|2) and (Q, determine a (strictly) feasible 
wake-up schedule. 

We denote as violation a condition such that constraints Q 
are not satisfied and yet a wake-up schedule is feasible. 

A. Key Assumptions 

We summarize here the main assumptions that we will use 
in the remainder of the paper: 

• Nodes do not have global knowledge of the common time 
axis; 

• Nodes do not have global knowledge of the Li and Ui 
values assigned to other nodes in the network; 

• Energy consumption is directly proportional to the duty 
cycle. This implies that we neglect non-linear effects and 
energy costs associated to wake-up operations; 

• Neighboring nodes are synchronized on common time 
slot boundaries ||2]- 

We remark that our results apply to general networks; however, 
due to the customary need to prolong battery lifetime, wireless 
sensor networks will be used as our reference case throughout 
the paper 

IV. Characterization of the Wake-Up Problem 

In this section we provide an algebraic characterization of 
the WAKE-UP problem. Using standard notation, we denote 
by (r, s) and by [r, s] the greatest common divider (gcd) and 
the least common multiplier (1cm) of r, s £ Z, respectiveljQ. 
Integer r divides s, which we write as r\s, if s = r ■ v for 
some integer v. 

We first recall a fundamental result: 



from which 



Theorem 1 ( Chinese Remainder theorem ^21'^ ) Let 
Til, 7i2, ai, 02 G 1i, consider the system 



X = ai mod ni 
X = a2 mod n2 



(3) 



Then, 

(i) system (O has a solution if and only if (rii , 712) |ai — a2; 
(iii) the general solution of ^ is given by x = Xq 
mod [^1,^2], where xq is a particular solution 0/(0. 

In order to derive a particular solution xq at step (ii), 
the extended Euclidean algorithm II2TI guarantees that integer 
Mi,M2 exist that solve the following Bezout's identity 

niui - 712U2 = 'T-2) 

' We use notation [r, s] to also denote the interval of the integers u such 
that r < u < s; the meaning of the notation is clear from the context. 



M i(ai - 02) 

("-l,"-2) 



02 - n2 - 



"2(01 



02) 



(»^l,"-2) 



The extended Euclidean algorithm has complexity 
0(log^ max{ni, 712}) = 0(1). It is clear that when 
system (|3]l consists of fc > 2 congruences, it can be solved 
iteratively; the solution can be determined repeating the 
operation for a single pair fc — 1 times. The time complexity 
for one step is 0(1), which in turn becomes on the order of 
0{k). 

From the model introduced in Sec. |lll] the solution of the 
WAKE -UP problem has to satisfy the following relation: 



lij + h ■ Ui ~ ttji + k 



(4) 



for ij € E and Oy G Ai and Oji G Aj, where h,k € Z. The 
sequence of rendezvous times {x2j}n£Z can therefore be seen 
as solution of ([3]i, with oi = o^j, rii = n,;, 02 = aji and 

Existence of feasible schedules follows as a consequence of 
the Chinese Remainder theorem. 

Theorem 2 A strictly feasible wake-up schedule always ex- 
ists, and can be constructed in time 0{N) with a fully dis- 
tributed asynchronous algorithm requiring 0{M) messages. 

Proof: A feasible wake up schedule is constructed in N 
steps choosing for all nodes rii > Li such that {ni,nj) ~ 1 
for all pairs neighboring nodes i and j. For each edge we 
can apply the Chinese Reminder Theorem to its end nodes. 
Irrespective of o,; and aj, it holds 1 = {ni,nj)\{ai — aj) 
at each step. The algorithm attains by construction a feasible 
wake-up schedule. ■ 

Remark 1 The above construction shows a basic effect: the 
request of synchronization of duty cycles introduces a delay 



drift (DD). In the worst case since [ni,nj] 



Tij, prod- 



ucts become exponentially large with respect to the network 
diameter. In order to minimize such drift, one would request 
[rii^rLj] = min{r7,i, rij} at each step, i.e., either ni \nj or vice 



We get therefore: 

Corollary 1 The Weak WAKE-UP problem can be solved in 
polynomial time. 

In the general case, when constraints on the maximum 
wake-up period are imposed, the problem becomes difficult 
due to the simultaneous requirements on the phases of nodes. 
Given a certificate wake-up schedule, we can verify with 2E 
congruences if the constraints (|2]i and ([T) are satisfied. Also, 
congruences appearing in (|4|i have the cost of computing 
{rii, nj) for each pair of nodes i.e., a solution of Strong 
WAKE -UP problem can be verified in polynomial time with 
respect to the input size, i.e., M. 

In the general case, when constraints on the maximum 
wake-up period are imposed, the given problem is made 



Algorithm 1 PER10D{L,,U,) 

Require: Li <Ui ^ Set : x = Li,ni = Li 
1: while X <Ui do 
2: if B\x then 
3: Hi ^ x; BREAK 
4: end if 
5: x <— a; + 1 
6: end while 



difficult due to the simultaneous requirement on the phases 
of nodes. It indeed belongs to the polynomial time verifiable 
algorithms class, i.e., WAKE-UP is NP. In fact, given a certifi- 
cate wake-up schedule, we can verify with 2E congruences 
if the constraints (|2]i and ^ are satisfied. Also, congruences 
appearing (|4|i have the cost of computing {ni,nj) for each 
pair of nodes i.e., a solution of WAKE-UP problem can 
be verified with 0{M) operations. 

The wake-up scheduling problem relates to the INTEGER 
FACTORIZATION problem: Given integer R, determine the 
factorization of R. 

Notice that while the answer to the question whether a 
factorization of R exists - namely the PRIMALITY prob- 
lem - can be actually solved in polynomial time Il22l . the 
INTEGER FACTORIZATION problem is harder At present 
there exist several algorithms for factoring an integer R 
that have estimated complexity either exponential in the 
size of the binary representation of R (as in the case of 
the exhaustive trial method), or sub-exponential (as in the 
case of the continuous fraction method 1*231). No polynomial 
time deterministic algorithm is known that solves INTEGER 
FACTORIZATION 1241 . 

Theorem 3 WAKE-UP is as hard as INTEGER 
FACTORIZATION. 

Proof: Consider a polynomial time reduction to the 
INTEGER FACTORIZATION problem. The input is a integer 
R > 0. The reduced instance of the WAKE -UP is a graph 
composed of two nodes, 1 and 2, and one edge that joins 
them; furthermore U2 = L2 = R, Ui = R — 1 and Li = 2. 
Assume that a solution of the INTEGER FACTORIZATION 
problem is given, i.e., R has proper factors m and n, i.e., 
R = m ■ n, m /I, nl /I. Consider n2 = L2 ^ U2 ~ R and 
ni ~ m. Furthermore, let ai = 02 = 0: then (0, m) and (0, R) 
are a strictly tight wake-up schedule. Vice versa, assume that 
a feasible schedule exists. Indeed n2 ~ L2 = U2 — R^y con- 
struction. Observe that the ni > 1 since Li > 2. However, the 
Chinese Reminder Theorem ensures that < R, which 

in turn implies that ni and R cannot be co-prime, otherwise 
[ni,R] = niR > 2R > R, a contradiction. Hence ni\R, and 
it is a proper factor of R since 1 < 2 = L2 < rti < R— 1 < R. 
Hence, by solving the WAKE-UP problem, we have answered 
to the problem of factorizing R as R — ni ■ q for some pair 
of proper factors ni,q > 1, which concludes the proof. ■ 



As reported in App. |D] there exists a further argument 
for the equivalence proved above: a connection exists be- 
tween WAKE-UP and INTEGER FACTORIZATION through 
a standard isomorphism. Also, WAKE -UP is not NP-complete 
unless NP = co-NP since it is known that INTEGER 
FACTORIZATION is NP-complete if and only if NP = co- 
NP, which (at present) is believed a convincing argument 
for ensuring that no polynomial time complexity algorithms 
exists Ii24\l . 

A question naturally arising is whether similar issues hold 
if we formulate the problem in continuous time, as opposed to 
the discrete-time framework. In App. |C]we show that, under a 
rationality assumptions, the two formulations are equivalent. 

Therefore, no efficient algorithm for exactly solving 
Strong WAKE-UP can be designed. In the next section 
we will introduce two heuristics for the Strong WAKE -UP 
which have message complexity linear in the number of links 
in the graph. 

A. Local time versus global time 

From the implementation standpoint, a remark is needed: 
according to the framework proposed, the elements of Ai are 
defined with respect to the common time. But, one might 
question that nodes operations are likely based on local time 
only. However, from the Chinese Remainder theorem, phase 
differences mod (ti;, Tij) matter In practice, nodes will 
know their relative phase, i.e., the time when a wake-up occurs 
for a neighbor with respect to their own time axis. Thus, node i 
will store the phase corresponding to node j as A^^ = — ai 
and node j will do the same for A^j = aji — aj. 

We can now see that phase differences with respect to the 
global time axis can be handled using values measured with 
respect to local time axes. In fact, 

ttij — aji mod {iii,nj) 

= {A^j - Ajj) mod {ni,nj) - (ai - aj) mod {ni,nj) 
= (Ay - Ajj) mod {ni,nj) + (a^ - ai) mod {ni,nj) 

Where last term corresponds to the offset of the two local 
clocks and can be easily measured time-stamping the local 
index of a certain time slot. 

In the rest of the paper we will refer to the common time 
axis for the sake of clarity. 

V. Algorithms Design 

We assume that each node i maintains a sorted list of 
neighbors Ng{i) according to their degree di. 

The algorithms that we adopt in the following use a certain 
factor basis B = {pi,p2, ■ ■ ■ ,Pn} where piS are prime 
numbers 1211 . For the sake of notation, let us say that B\n 
if n factorizes in B. Observe that the choice of the factors 
of the wake-up periods at each node is critical; the Chinese 
Remainder Theorem suggests that period rii should be chosen 
free of large factors. 

Nodes choose period rii € [L,, Ui] such that it factorizes 
according to B using the procedure described in Alg. [T] Alg. [T] 



Algorithm 2 BFS WAKE-UP 

Require: Li <Ui ^ Default : Ai = {ai} 
1: for all i e F do 

2: rij = PERIOD(Lj, Ut) 
3: fit) ^ 

4: end for 

Ensure: Select: i ^ V with largest degree 

5: fit) = 1 

6; Q {i} {Enqueue the first node} 

7: while Q 7^ {Visit all the nodes in BFS} do 

8: i <r- Q[l] {Get the node on top of the queue} 

9: for all j e Ng{i) do 

10: if f{j) ~ {If it was not explored before} then 

11: Q <r- {j} {Enqueue this node} 

12: aj ■(— ai {Align time axis of node i and node j} 

13: Aj <r- {uj} {Update the phase if j} 

14: Xij = Xji ^ ai + k-[ni,nj], A: G Z {Rendezvous 
times} 

15: fij) ^ 1 {Visited} 

16: end if 

17: end for 

18: Q ~ Q \ {i} {Dequeue i} 

19: Tii ^ [rii, (riij , . . . , Ui^ )] {Use largest possible period 
at i} 

20: f{i) = 2 {Explored} 

21: end while 



returns the smallest integer in [Li, Ui] that can be written as 
pV ■ ■ 'Pn"' < e Z, provided that such a number exists. 
If it does not exist, Li is returned. We now introduce two 
heuristics for Strong WAKE-UP. 

A. BFS WAKE -UP Algorithm 

Once nodes have chosen their periods according to Alg. [T] 
the simplest rationale is to align their phases pairwise. Alg. |2] 
is based on a breadth first search (BFS) visit of the network 
rooted at a given node; as customary for BFS search a queue 
Q is employed. 

Each node i maintains information about nodes in its 
neighborhood, i.e., the node ID, the phases set Ai and a special 
flag / that is used for the BSE procedure, where (i) f{i) — 
means j unexplored (ii) /(j) = 1 means node j discovered 
but not all neighbors of j examined (iii) /(j) = 2 node j 
explored. By default, node i wakes up at the origin of its own 
period, i.e., at a,;. 

The algorithm selects all unexplored nodes and enqueues 
them at most once, so that it terminates in N steps. 

Alg. |2] reports on the pseudocode of the centralized ver- 
sion of the BFS WAKE -UP algorithm; however, the BFS 
WAKE -UP algorithm can be easily implemented in a dis- 
tributed fashion [251. Observe that BFS WAKE -UP attains 
phase synchronization and a strictly feasible schedule. 

Below, we formalize the properties of the algorithm; the 
proof is reported in App. lAl 



Theorem 4 BFS WAKE -UP has the following properties: 

i. it produces a strictly feasible wake-up schedule; 

ii. it has complexity 0{E); 

iii. a distributed implementation requires 0{M) messages. 

Remark 2 For implementation purposes, the above algorithm 
can be very easily coupled to the construction of shortest 
path trees i l25l/ ; in such a way customary algorithms such as 
Dijkstra can conveniently provide joint routing and wake-up 
scheduling with no need for separate message exchange; in 
turn, the weight of link ij is represented naturally by [n^, nj]. 

The convergence time of the algorithm is proportional to the 
diameter of the network |25]; the BFS WAKE -UP requires a 
preliminary leader election procedure; further its convergence 
may be slow. Accordingly, we devised an alternative heuristic 
algorithm, described below, which effectively addresses such 
shortcomings by parallelising the computation of the schedul- 
ing. 

B. PARALLEL WAKE-UP Algorithm 

The design of a parallel algorithm for wake-up synchroniza- 
tion requires careful handling of nodes phases; in particular, 
in order to avoid the increase of the size of the AiS, it is 
possible to operate pairwise adaptation of nodes phases (phase 
adaptation). 

Lemma 1 Let ij,jk £ E and fix ni,nj,nk, auj and aij. A 
feasible wake-up schedule for i,j,k such that ajk = dji exists 
if and only if (rij, (ni,nk))\{akj —aij). If such condition holds, 
the schedule is obtained for Ujk = aji = akj + n-k • (xq ■ e), 
where e,Xo G Z are solutions of akj — aij = (rij , Uk) • e ■ Xq -\- 
{ni,nj) -e-yo, with ya e Z. 

The proof is reported in App. |A] 

The algorithm pseudo-code is reported in Alg. |3] Each node 
i maintains information about nodes in its neighborhood, i.e., 
the node ID, the phases set Aj, the wake-up period rij and 
two special flags / and t. Special flags are used in order to 
determine Ai and n^. In particular (i) f{j) ~ means j 
unexplored (ii) f{j) = 1 means node j discovered but rij has 
not been fixed yet and all neighbors of j have been explored 
(iii) f — 2 when node j has been explored and therefore aj 
and Tij have been fixed. Also, t{j) = means j that can adjust 
aj (ii) t{j) = 1 means i cannot change aj value because some 
of its neighbors used aj. 

At the start of the algorithm each node i chooses ri; using 
PERIOD and maintains a list of neighbors Ng{i). Then node i 
is marked unexplored and available to change ai value {f[i] ~ 
t[i] = 0). Hence, each node starts a random timer: when the 
timer expires node i freezes neighbors' timers and determines 
rendezvous times with all nodes in its neighborhood. 

If i is unexplored, i.e., f{i) = 0, node i checks whether it 
is in range of an already existing node j explored. If so, i just 
aligns ai = aj; j is then is marked not available to change its 
time axis anymore, i.e., t{j) = 1. Otherwise, if i is discovered, 
i.e., /(j) = 1, i checks whether it is in range of two already 



Algorithm 3 PARALLEL WAKE-UP 

Require: Li < Ui 
1: for all i G F do 

2: rii = PERIOD(Li, Ui) Default : A; = {aj 
3: f{i),t{i) <— {Fix the initial value to flags / and t} 
4: end for 
Ensure: Node i starts a random timer 
5: Timer expired at node i: send a message to stop neighbors' 

timers in Ng{i) 
6: Initialize: 

7: Tii ^ [ui, (riii J • ■ • > "-id. )] {Use the largest possible period 
at i} 

8: if f{i) = and there exists j e iVg()i s.t. f{j) ^ 
then 

9: a; 4- Q!j {Assign the neighbour's phase to node i} 
10: t(j) ^ 1 {Node j cannot change its phase anymore} 
11: else if there exist j,k e Ng{)i : f{j) ~ f{k) = 2 and 
Aj n Afc = and conditions of Lemma [T] are satisfied 
then 

12: case tii) = then ^ ajk 

13: case t{i) = 1 then Ai = A^ U {aju} 

14: A; \ {aij.aik} {Delete and aiu from AJ 

15: end if 

16: Visit all neighbors: 

17: for all i e Ng{i) and f{j) ^ 2 do 

18: if /(j) = {If it was not explored before} then 

19: aj ai {Assign the phase to node j} 

20: f{j) = 1 

21: else if At n A^- = then 

22: Aj Aj U {ai} {Add the phase to node j} 
23: end if 

24: x^j = ai + u {u computed using extended 

Euclidean algorithm} 
25: Xij — Xji = x^j + rij], fc = 0, 1, . . . {Rendezvous 

times} 

26: /(j) = 2 

27: end for 



existing j and fc already explored without a phase in common 
between them and if the phase can be adjusted then i has two 
options: (i) if i is available to change its phase {t{i) = 0), 
then i just aligns its phase to the novel phase, i.e., ai = aj, 
(ii) if i is not available to change its phase {t{i) = 1), then 
the novel phase is added to Ai. 

Then, i visits its neighbors for which no adaptation was 
performed. If j g Ng{i) is unexplored, aj = ai and j is 
marked discovered, f{j) = 1. If j is explored and AiflAj = 0, 
the phase of i is added to Aj. Finally, i marks itself explored. 

The proof of correctness and the fact that Alg. |3] has 
message exchange is 0{M) is similar to the one for Alg.|2l 

Theorem 5 PARALLEL WAKE-UP.- 
i. produces a wake-up schedule; 
iii. it requires 0{M) messages for a distributed implemen- 
tation. 



In principle, since multiple phases may be added at one 
node, it is possible that the output of the algorithm is a 
schedule that is not feasible. A simple bound on the average 
duty cycle (DC) is provided by the following 

Theorem 6 Under Alg.\3\DC < where L = rainLi and 
d is the average node degree 

Proof: The algorithm adds at most di — 1 phases per node; 

Remark 3 Observe that the above algorithm can be used to 
adapt the wake-up scheduling dynamically in case one or more 
nodes join the network. In particular, it is possible to combine 
the two algorithms in order to initialize a wake-up sched- 
ule network-wide with BFS WAKE-UP, and use PARALLEL 
WAKE -UP to add new nodes to existing schedules or to repair 
locally a schedule, should a node fail or loose synchronization. 

C. Example 

We reported in Fig. |2] a simple example illustrating the 
BFS WAKE-UP and of the Parallel WAKE-UP algorithm 
on a sample graph with 7 nodes. In that instance U = 20, 
L, = (2,3,9,7,11,5,2), a, = (1,3,7,9,6,1,0) and B = 
{2}. Several wake-up periods coexist: 4, 8, 16. 

The BFS WAKE-UP run generates the list of visited nodes 
(1, 2, 3, 4, 5, 6, 7); observe that every node aligned to ai = 1. 
In the Parallel WAKE-UP node 2 and 7 perform the 
algorithm operations first, then 1, 5, 6 and finally 4. As we 
can see from the example. Parallel WAKE-UP assigns 
heterogeneous phases; at node 4 phase adaptation according 
to Thm. [T]is not possible so that finally A4 = {1,3}. 

D. Particular cases 

There exists specific cases when strictly tight wake-up 
schedules are attained easily. 

Lemma 2 Let for any node i be B\Li, B\U = Ui and Li < 
U, then there exists a strictly tight wake-up schedule; it can 
be attained using BFS WAKE-UP. 

Proof: Consider any link ij and let rii — Li and 
rij = Lj . The statement follows immediately from the Chinese 
Reminder theorem since [ni,nj]\U and one can use BFS 
WAKE UP to assign the same phase to all nodes in the network. 

■ 

In particular the following example shows a simple applica- 
tion of the above case. Consider a tree topology T and assume 
nodes are labeled in order to respect the partial order induced 
by the distance from the root node 1. Let Lo be associated 
to the root node, and Lk = L^p^ be associated to nodes at 
level k. A feasible schedule can be constructed as follows; 
root node wakes at times xq ~ mod Lf). Node j wakes 
up at Xj ~ ao mod p^^'^^^Lq where rj is the level occupied 
by node j and g is any positive integer function. In this case 

U = LQp"'''''9(r,)^ 
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Fig. 2: Example. BFS WAKE-UP and PARALLEL WAKE-UP on a sample topology. The table reports on the final output of the algorithms. 



VI. Performance Evaluation 

In this section, we report on the outcomes of experiments for 
the BFS WAKE-UP and PARALLEL WAKE-UP algorithms. 
In the first set of experiments, performed using Matlab®, we 
tested the algorithms on random topologies. 

We considered a square playground of side 100 m with 200 
nodes randomly deployed according to a uniform distribution. 
Nodes were assumed to be connected if their mutual distance 
was below a given value (communication range), varied in the 
range 10 — 80m. For each setting, 100 runs were performed. 
For the construction of the schedules, two factor bases were 
considered, B = {2} and B = {2,3,5}. The constraints were 
set by taking, for each node, Li uniformly distributed in [1 , 35] 
and Ui uniformly distributed in [50, 100]. 

We considered as performance metric the average duty 
cycle, as introduced in the previous section, and the average 
normalised delay drift, defined as: 

DD = ^y y (5) 

It is worth remarking that DD > 1 implies that some 
constraints on the delay were violated. 

Results are reported in Fig. |3] In terms of duty cycle, BFS 
WAKE-UP outperforms PARALLEL WAKE-UP, as expected, 
since the second one may add multiple phases per node. The 
performance difference reduces in dense networks (i.e., with 
large transmission range). A closer look at the results showed 
that, for the parameters chosen, the schedules generated by 
PARALLEL WAKE-UP were not feasible in less than 1% of 
the cases. In terms of average normalised delay drift, the 
two algorithms attain the same performance. BFS WAKE -UP 
is much slower in converging in sparse networks, but its 
performance increase (and, actually, outperforms PARALLEL 
WAKE-UP ) in very dense networks. As it may be seen, the 
use of a richer basis improves significantly the performance 
in terms of duty cycle, but it also leads to worse delay drift 
performance. 

In general, we may conclude that PARALLEL WAKE-UP 
represents a meaningful choice in most situations, presenting 
performance close to BFS but with a shorter convergence time. 



Violation 


B 


[75, lOOJ 


[60, 100] 


[45, 100] 


Magnitude 


{2| 





2.04 ±0.07 


9.32 ± 0.26 


Magnitude 


{2,3,5} 








88.38 ± 1.84 


Percentage 


|2| 





1.25% ± 0.12 


5.43% ± 0.27 


Percentage 


{2,3,5} 








9.36% ± 0.43 



TABLE II: Perfoimance of BFS WAKE-UP and PARALLEL 
WAKE-UP, fraction and magnitude of violations. Settings as in Fig.|3] 



Discrete events simulation: In a second set of experiments, 
we implemented the proposed techniques in an event-based 
simulator, Omnet++, to assess the advantage of supporting 
heterogeneous wake-up schedules. We considered 50 sensors 
deployed over a IOO7T1 x lOOm area, with communication range 
of lOm and set the time slot duration to lOOms. We assumed 
that packets were generated at each node every 30 s, and that 
all packets were routed towards a sink, located in (1, 1). A tree 
was constructed at the beginning of the simulation to route 
packets to the sink. 

As sensors far from the sink relay less traffic, their bounds 
on energy consumption can be tightened. We considered a 
situation in which the constraints Li depended on the distance 
from the sink as follows: Li = 2 for nodes at distance < 2 
hops, Li = 4 for nodes at distance of 3 and 4 hops and 
so on. As performance metrics we considered the average 
packet delay (from source to sink) and the lifetime of the 
nodes. The latency was computed up to the time at which 
the first node died. We compared the performance attained by 
BFS (PARALLEL turned out to offer the same performance) 
with those obtained in the absence of sleep mode ('No Power 
Saving') and those obtained when all nodes go to sleep 
synchronously for one slot every two ('Uniform'). The results 
are reported, plotted against the distance from the sink, in 
Fig. |4] As it can be seen, BFS is able to effectively prolong 
the lifetime of sensors far from the sink, while at the same 
time limiting the increase in packet latency. 

VII. Conclusions 

In this paper we introduced a framework for the synchro- 
nization of wake-up schedules, which plays a central role for 
energy saving in wireless networks, under the joint request 
to satisfy per node energy and delay constraints. We showed 
that the related WAKE -UP problem has a strong relation with 
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Fig. 4: Discrete event simulation: a) average node lifetime as a 
function of the distance from the sink, b) average delay. 
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Fig. 3; Perfomance of BFS WAKE-UP and PARALLEL WAKE-UP: 
results for Li G [1,35] and Ui G [50, 100]: a) average duty cycle 
(DC) b) average delay drift (DD) and c) convergence time. 



integer factorization, which adds a novel algebraic perspective. 

This work left out several interesting directions. In particu- 
lar, for relatively small numbers, factorization can be done in 
practice ||23l|. In turn, the optimal choice of the factor basis 
and of periods is a key problem in order to limit the drift of 
the rendezvous times. 

Furthermore, mapping end to end delay onto bounds on 
the wake-up period has been only touched briefly in App. |E] 
However, while there exist several works in literature tackling 
the general problem of bounding the end-to-end delay ll26l . 
1271 . a network calculus accounting for power saving duty 
cycles has not been proposed so far. Relating the WAKE -UP 
problem to those works is part of future work. 
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Appendix 

A. Proof of Thin. ^ 

Proof: i. Consider node i and assume that the algorithm 
started at node 1 for simplicity: at the end of the algorithm 
it holds {Ai,n.i) = ({aj, [(n^, n.^ J, . . . , (n^, Notice 
that tti — aj- ~ for all ji G Ni and also 

([n,,njj, . . . , = [n^, (n^,, . . . , (6) 

Hence, every pair {ai,ni) solves the modular congruence 
system according to the Chinese Remainder theorem. 

ii. For each node j, Ng{j) operations are performed, 
including the calculation of [nuj] , which requires 
on 0(log'^ max(rii, nj)) operations. The overall 
complexity of the algorithm is then bounded as 
E»=i J2JeN,{^) log"^ max(n„ n^) = 2E\og^ U 

iii. The statement follows immediately by a double counting 
argument since every node i needs to exchange at most 2 
messages with each neighbor in order to visit neighbors, 
communicate its phase and receive acknowledgment. ■ 

B. Proof of Lemma Q] 

Proof: A novel phase a can be given to node j iff the 
following system solves 

X — aij mod rii ^ a mod nj 

= akj mod Uk (7) 

from which {nk,nj)\akj — a and {ni,nj)\a — aij. Thus, if a 
novel phase can be assigned then a solution exists for 

c = - a^j = X {nj,nk) + y {ni^rij), x,y 

From O, d ~ (ni,{nj,nk))\akj — ciij'. let de = a^-j — aij 
consider a solution , j/o 

d= {nj,nk)xo + (ni,nj)yo, xo,yoeZ 

then 

de^ {nj,nk)exo + {n^,nj)eyo 
the novel phase writes 

a = akj - {nj,nk)exo 

Observe that a — aij = akj — {rij, nk)exo — aij ~ {rii, nj)eyo 
which concludes the proof. ■ 

C. Continuous-Time Analogue 

Here we consider wake-up scheduling in continuous time, 
under the additional requirement that the difference of phase 
at different nodes is rational (observe that this is a practical 
requirement since numerical representation is rational). We 
conclude that the continuous time approach plus the request 
for a rational difference of phase at two nodes is equivalent 
to operate wake-up scheduling in the discrete time domain. In 
fact, consider the general case of a network where wake-up 
scheduling is operated and nodes are assigned phases ai and 
period Ti. It holds for any node i: 

KijTi + = KjiTj + aj, j G Ng{i) 



Now, assume G is connected, and consider TV paths 
Pi, P2, ■ • ■ 7 Pn from node 1 to nodes j ~ 1, . . . ,N, respec- 
tively: it is easy to see that the following set of equalities 
hold 

Tj/Ti = qji + ^ PrsClrs 
rsePj 

where a^s = ctr ~ cts, whereas qji and prs are rational. Thus, 
it follows that also Tj/Ti is rational. This is equivalent to 
admit that the overall system time, i.e., the period of each 
node Ti can be reduced to a convenient multiple of a certain 

To = Ti/TVo. 

D. Standard isomorphism 

A basic result in algebra that provides a natural link of the 
Chinese Remainder theorem with the problem of factorization 
of integers. With standard notation denote Z/nZ the quotient 
ring of integers modulus n 1281 . Define the standard isomor- 
phism 

: Z/nZ X Z/mZ — > 'L/[m, n]Z 
(a, 6) -> c 

where c solves for 

c = a mod n 
c = b mod m 

Assume now that m, n are two coprime integers, i.e., (to, n) = 
1; it is easy to see that ^ is an isomorphism. Translating with 
the terminology used in wake-up scheduling used through- 
out this paper, this means that for every wake-up schedule 
{(a, n), {b, m)}, there exists a unique rendezvous time cq such 
that c = Co mod [to, n] for any other rendezvous time c 
generated by the schedule. Viceversa, given a non-empty set of 
rendezvous times for two nodes having coprime periods, there 
exists a unique pair of phases a and h that corresponds to such 
set. This also means that phases are automatically determined 
by the set of rendezvous times and the coprime periods, 
i.e., the only information needed is that mn = [m, n], i.e., 
(m, n) = 1. This provides another interpretation of the reason 
why WAKE-UP is as hard as INTEGER FACTORIZATION, 
as proved in Sec. ITV] 

E. On the Dimensioning of the Upper Bounds 

In general, application-level constraints for the application 
setting we target are expected to be given in terms of: 

> Expected lifetime of the network. This can be translated 
in terms of an upper bound on the duty cycle of single 
nodes. This is reflected in a constraint of the form like 

> End-to-end delay. Bounds will be given in terms of the 
delay incurred between the generation of a message at 
a node (which corresponds, in our framework, to the 
reading of a sensor) and the delivery at the appropriate 
sink, where such data may be processed or transferred, 
via some other communications technology, to a remote 
processing center 



For the latter case, we need to identify means for translating 
end-to-end delay constraints into node-level delay constraints 
of the form (|2|i. The problem is then to translate a number of 
constraints, given on routes, to constraints on the single links 
constituting the routes. 

In order to do so, we start by introducing some notation. 
We define by Ri the route starting from node i and ending 
in one of the sink nodes. We say that node j lies along the 
route Ri if all traffic originating from i passes through j in 
order to reach the sink, and write j E Ri. The length (in hops) 
of route R{i) is denoted by \R{i)\- Conversely, we denote by 
the set of routes passing through node i. The worst-case 
delay associated to route Ri is denoted by D{Ri). The delay 
constraint on route Ri is expressed as D{Ri) < rji, where rji 
is a constant. 

The difficulty of the problem lies in the fact that the 
constraint on the energy expenditure dl) impacts the minimal 
latency experienced by a message passing through a given 
node. 

We assume that (i) routes are static throughout the lifetime 
of the network (ii) nodes arriving during an active period are 
sent at the next active period occurrence (iii) the duration of 
a time slot is sufficient for transmitting all messages queued 
at any given node. 

Now we consider the latency experienced at a given node 
and imposed by the duty cycle. Under the aforementioned 
assumptions, the worst-case delay incurred by a message 
passing through node i can be lower bounded by + 1. 
The term 'worst-case' here refers to two factors. One, to the 
arrival time of the message, whereby the worst case coincides 
with a message arriving right after the beginning of an active 
slot. Second, to the rendez-vous points. In order to be able to 
move the message one hop further along the route, it is not 
sufficient for node i to be in the active state, but, given the 
system model presented in Sec. |III] it requires also the next 
hop node to be in the active state. In general, the delay will be 
a multiple of plus the time slot necessary to forward the 
message, iven condition ([T]), this implies that such worst- 
case delay is lower bounded by Lj + 1. 

Clearly, a set of feasible Ui has to satisfy Ui > Li + 1 for 
all nodes i. Furthermore, in order to be compatible with the 
end-to-end delay constraints it has to be J2jeB — ^'^^ 
all routes i?, . 



^It is worth remarking that the bound is exact if \Ai \ = 1. For \Ai\ > 1 it 
is an approximation of the actual one, derived assuming that active slots are 
uniformly distributed within the wake-up period. 



We thus need to find a solution of the following system of 
inequalities: 



jeRi 



jeRN 
Ui > 1 + 



(8) 



Un>1 + L 



The system has a solution if and only if the following 
condition is satisfied; 



,iV. 



(9) 



If such a condition is respected, a feasible solution can be 
found as follows. We associate to route Ri a "delay budget" 
rji. We consider how much is consumed by the duty cycle at 
the various nodes along the route Ri, which, using the bound 
outlined above, will be taken as + Ejej? -^j- '^^en 
distribute the remaining "delay budget" in a uniform manner 
across all nodes of the route Ri. We do this for all routes. Then, 
for each node i, we consider as Ui the minimum value of such 
quantity among all routes passing through it. In symbols: 



Ui ~ Li -\- 1 + mill 



\LL\ 



(10) 



It is trivial to prove that under (|9) solutions of the form 
([Tol l satisfy Ui > ii + 1 Vi. In order to prove that the other 
conditions in ([8]) are satisfied we note that: 



EUi= < + 1 + mm 



m-\Rk\- E Lh 
\Rk\ 



> < 



ieRj ieRj 



(11) 

Vj - iRjl - E Lh 

h^Rj 



(12) 

= Y,L^ + \Rj\ + iij - \Rj\ - ^ L,, = 77j, 

ie-Rj h£Rj 

(13) 

where we used the fact that if i G Rj then Rj G Fi and that 
minjxi} < Xj Vj G S. 

We recall that the method outlined above provides only 
one of the possible solutions of the system dHJ. In particular, 
there is no reason for which a uniform distribution of the 
remaining delay budget along a route is optimal, though this 
seems reasonable in a variety of settings. 

For the special case in which Li ^ LMi and rji = 77 Vi, i.e., 
all nodes have the same energy constraint and all messages 



should be delivered to a sink within the same deadline, 
regardless of their distance, ( \T0\ simplifies as: 



Ui = L + 1 + mill 



'7- l^il -E/ceH, L 



= L + 1 + min 



V~\R,\-{L + l) 



V 



max iR-j 



Rj eJ^i 



(14) 

(15) 
(16) 



In the uniform case, therefore, the bounds Ui depend only 
on the maximal length of the routes passing through node i. 



